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The Skyrme-Faddeev-Niemi (SFN) model which is an 0(3) a model in three dimensional 
space up to fourth-order in the first derivative is regarded as a low-energy effective theory 
of SU(2) Yang-Mills theory. One can show from the Wilsonian renormalization group 
argument that the effective action of Yang-Mills theory recovers the SFN in the infrared 
region. However, the theory contains another fourth-order term which destabilizes the 
soliton solution. In this paper we derive an extended action including second derivative 
terms and obtain soliton solutions numerically. A new topological lower bound formula 
is infered for the extended action. 
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1. Introduction 

The Skyrme-Faddeev-Niemi (SFN) model which is an 0(3) cr model in three dimen- 
sional space up to fourth-order in the first derivative possesses topological soliton 
solutions with torus or knot-like structure. The model was initiated in 70's^and its 
interest has been extensively growing. The numerical simulations were performed 
l ^h/ l ^ l v l O I^ the integrability was showrt^, and the application to the condensed mat- 
ter physics ^ and the Weinberg-Salam model ^ were also considered. The recent 
research especially foc uses on the consistency between the SFN and fundamental 
theories such as QCD In those references, it is claimed that the SFN 

action should be derived from the SU(2) Yang-Mills (YM) action at low energies. 
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One can also show from the Wilsonian renormaUzation group argument that the 
effective action of Yang-Mihs theory recovers the SFN in the infrared region li^. 
However, the derivative expansion for slowly varying fields n up to quartic order 
produces an additional fourth-order term in the SFN model, resulting in instability 
of the soliton solution. 

Similar situations can be seen also in various topological soliton models. In the 
Skyrme model, the chirally invariant lagrangian with quarks exhibits fourth order 
terms after the derivative expansion and they destabilize the soliton solution^^l 
To recover the stability of the skyrmion, the author of Ref J7. introduced a large 
number of higher order terms in the first derivative whose coefficients were deter- 
mined from the coefficients of the Skyrme model by using the recursion relations. 
Alternatively, in Ref . El Gies pointed out the possibility that the second derivative 
order term can work as a stabilizer for the soliton. Similar form o f th e action was 
also proposed by Forkel with somewhat different expansion schem^^. The author 
also found out the saddle point soliton solutions with hedgehog type. This result 
is quite encouraging to examine the consistency between the topological soliton 
physics and the Yang-Mills theory. 

In this paper, we compute the extended Hopf soliton solutions from the action 
proposed by Gies. In section El we give an introduction to the Skyrme-Faddeev- 
Niemi model with its topological property. In section |3 we show how to derive the 
SFN model action from the SU(2) Yang-Mills theory. In sectional soliton solutions 
of this truncated YM action are studied. For this purpose, we introduce a second 
derivative term which can be derived in a perturbative manner. The naive extrem- 
ization scheme, however, produces the fourth order differential equation and the 
model has no stable soliton solution. Failure of finding the soliton is caused by the 
basic feature of the second derivative field theory. In section|31 the higher derivative 
theory and Ostrogradski's formulation are reviewed. We show the absence of the 
energy bound in the second derivative theory using an example in quantum me- 
chanics and introduce the perturbative treatment for the second derivative theory. 
In section |H| we present our numerical results. The possibility of finding new topo- 
logical bound for this extended, perturbative soliton solutions is also discussed. In 
section [71 are concluding remarks. 



2. Skyrme-Faddeev-Niemi model 

The Faddeev-Niemi conjecture for the low-energy model of SU(2) Yang-Mills theory 
is expressed by the following effective action: 

5sFN d^x\^{d^nf + |-(9^n x d.nf] (1) 

where n{x) is a three component vector field normalized as n • n = 1. The mass 
scale A is a free parameter and in this paper we set A = 1. It has been shown that 
stable soliton solutions exist when gi > 0. 
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The static field n{x) maps n : i-^ and the configurations are classified by 
the topological maps characterized by a topological invariant H called Hopf charge 

H^^jA^F^F^dA (2) 

where F is the field strength and expressed in terms of n{x) as F = {n ■ dn A dn). 
The static energy Estt from the action Q has a topological lower bound , 

Estt > KH""^ (3) 

where K = Wir^y/gi. 

Performing numerical simulations, one can find that the static configurations for 
H = 1,2 have axial symmetry S Thus "the toroidal ansatz" is suitable to impose 
on these configurationl^. The ansatz is given by 



ni = v^l — /?) cos(7Vq! + v{r], (3)) , 

^2 = ^1 -w'^iViP) sin{Na + v{ri,l3)) , (4) 

where {rj, /?, a) is toroidal coordinates which are related to the Gaussian coordinates 
in as follows: 

a sinh 77 cos a a sinh 77 sin a a sin /? , , 

x = ,y= ,2 = (5) 

T T T 

with T = cosh 77 — cos/3. The function w(ri, (3) is subject to the boundary conditions 
w(0,/3) = l,w(oo,/3) = —1 and is periodic in (3. v(ri,l3) is set to be v{rj,f3) = 
Mf3 + vo{ri, (3) and uo(, (3) is considered as a constant map. Equation ^ then gives 
H = NM. 

To obtain soliton solutions with higher derivative terms, we propose a simpler 
ansatz than Q), which is defined by 



ni = ^Jl -w^{r]) cos{Na + M/3) , 

712 = ^/^-w'^ir]) sm{Na + MP) , (6) 

713 = w{t]) , 

where w(j]) satisfies the boundary conditions w(0) = 1,k;(oo) = —1. By using lO, 
the static energy is written in terms of the function ^(77) as 



Estt =2tt a dr] 



+ (1 - w^)Um.n{v) + ^ sinh7;cosh77(w') C^M,iv(??) 



'j' Um.n{v)= {M^ + - 

dn V SI 



drj' \ sinh^77 

The Euler-lagrange equation of motion is then derived as 

" /2 

" + UM,Niv)w + ^ (-2N^ coth^ Tjw' 
1^ 2a'^ \ 



WW 



1-w^ {1-w'^y 
+ (cosh^ 7/ + sinh^ 'n)U m ,n {11)^' + siTLh.'qcos\irjUM,N{'rj)w''^ = . (7) 
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The variation with respect to a produces the equation for a. Sohton solutions are 
obtained by solving the equations for a as well as for w. 

We obtained soliton solutions numerically for both ansatz Q and ©. The total 
energies together with the topological lower bound Q are shown as a function of 
coupling constant gi in Fig^ We found that this simple ansatz produces at most 
10% errors and it does not affect to the property of the soliton solution. 



3. Cho-Faddeev-Niemi-Shabanov decomposition and the effective 
action in the Yang-Mills theory 



In this section, we briefly review how to derive the SF N effec tive action from the 
action of SU(2) Yang- Mills theory in the infrared limit ll^lSl F or the gauge fields 
A^, the Cho-Faddeev-Niemi-Shabanov decomposition is applied ^^^^2lEI 



= nC^, + (9^n) xn + W^,. (8) 

The first two terms are the "electric" and "magnetic" Abelian connection, and W ^ 
are chosen so as to be orthogonal to n, i.e. W ^ • n = 0. Obviously, the degrees of 
freedom on the left- and right-hand side of Eq.JHl) do not match. While the LHS 
describes Scoior x 4Lorcntz = 12, the RHS consists of (C^ :)4Lorcntz + (ri :)2coior+(W^/i ■ 
)3coior X 4Lorentz — '^n W ^^=0 — 14 dcgrccs freedom. Shabanov introduced in his paper 
^1 the following constraint 

C^, W^) ^ 0, with X • n = . (9) 
The generating functional of YM theory can be written as 

Z = j VnVCVW5{x)^FP^se''^™~^<'' . (10) 

App and S^f are the Faddeev-Popov determinant and the gauge fixing action re- 
spectively, and Shabanov introduced another determinant As corresponding to the 
condition X = 0. YM and the gauge fixing action is given by 

/r 1 1 1 

Inserting Eq.||Sl) into the action, one obtains the vacuum functional 



with 
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Fig. 1. The total energy of the sohtons with (a) H = 1, (b) H = 2, as a function of the couphng 
constant gi: the simple ansatz JHJ, the toroidal ansatz J3] by Gladikowski-Hellmund, and the 
expected topological lower bound l3l . 



2, 2008 4:36 WSPC/INSTRUCTION FILE sawado 



6 N. Sawado 



and 



K 



d^{n X d^n) , (in gauge ctg = 1) 



The classical action for n including the gauge fixing term is given by 



Sc.] — 



I 



I 



{d^n X nf 



The S functional is expressed by its Fourier transform 

Integrating over "fast" variables C, W, <p, we finally obtain 

g-Se« ^ e-^-AFpAs(detM^)-i/2(detM^)-i/2(det-g^(M^)-ig^)-i 



(13) 
(14) 

(15) 

/2 

(16) 



where := i{—dfj_+d^n®n). Here several nonlocal terms and the higher derivative 
components have been neglected. 

Performing the derivative expansion for the res ultant determinants with dn, one 
obtain the effective action for "slow" variable 

rl 



{d,nf + ^(9^n X d,nf + f (a^n)^ 



92 , 



(17) 



For gi > and g2 = 0, the action is identical to the SFN effective action 

In order for soliton solutions to exist, g2 must be positive El. However, 52 is 
found to be negative according to the above analysis. Therefore we consider higher- 
derivative terms and investigate if the model with the higher-derivatives possess 
soliton solutions. 

4. Search for the soliton solutions (1) 



The static energy is derived from Ea. pTfl as 



'^(5.n)2 + |(5.nx9,n)2 + |(5.n)4 
E2{n) + E^l\n) + Ef\n) . 



(18) 



A spatial scaling behaviour of the static energy, so called Derr ick's scaling argument, 
can be applied to examine the stability of the soliton EOl. Considering the map 
X x' — fix (/i > 0), with n^'^' = n{^x), the static energy scales as 

e(Ai) = i;stt(n(^)) 



-E2{n) + ^i{E'i\n)+E\^>{n))^ 



p(2)/ 



(19) 
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Derrick's theorem states that if the function e{fi) has no stationary point, the theory 
has no static solutions of the field equation with finite density other than the vac- 
uum. Conversely, if e(/i) has stationary point, the possibility of having finite energy 



a) 



E 



(2)^ 



soliton solutions is not excluded. Ea. H19|) is stationary at /i = y E2/{E, 
Then, the following inequality 

gi{din X dju)"^ + g2(9in)^(5jn)^ 

^ 52(a,n • dju)^ (•.■ (9,n)2(a,n)' ^ (dm ■ d,nf) 

ensures the possibility of existence of the stable soliton solutions for g2 ^ 0. As 
mentioned in the section .92 should be negative at least within our derivative 
expansion analysis of YM theory. 

A promising idea to tackle the problem was suggested by Gies ' ' . The author 
considered the following type of effective action, accompanying an second derivative 
term 



9i 



92 , 



92, 



-{d^nY + '-^{d^n X d^nY - '-^{d^nf + • d'n) .(20) 

z o o o J 

Here we choose positive value of g2 and assign the explicit negative sign to the 
third term. In principle, it is possible to estimate the second derivative term by the 
derivative expansion without neglecting throughout the calculation. The calculation 
is, however, very laborious and hence we will show in detail somewhere else. 
The static energy of Eq. (|20|) with the ansatz ^ is written as 



Es\ 



2tt a j drj 



K) 



l\2 



{1-w'')Um,n{v) 



92 



sinh rf cosh 77 



'^2 



.91 
2a2 

2 



sinh ?7 cosh ?7(w') Um,n{^) 



+ {l-w')UM,N{ri) 



K) 



+ ^coth T] + sinh^ rj — sinh // cosh rj 

— 



w{w') 



'\3 



sinh 77 cosh 



11 — I 



{w'f 



1 — 

2 



v'Um,n{vi)^ 



+ ww" + {1-w'^)Um^n{v)\ +{w"f] 



where w" ^ 



The Euler-Lagrange equation of motion is derived by 



dE,, 



+ 



d fdEs, 



drf' V dw" / drj\ dw' 



dE^ 
dw 



= 0, 



(21) 



(22) 



which is too complicated to write explicitly. Thus we adopt the following notation 
fo{w,w',w") + gifi{w,w',w") g2f2{w,w' ,w" ,w^^\w^^'') = 0. (23) 
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Here w^^\w'^'^'^ represent the third and the fourth derivative with respect to rj. The 
first two terms of Eq. (|23|1 arc identical to those in Eq. Q . 

In addition to the boundary conditions at the origin and the infinity lo{Q) = 
I,a;(oo) = —I, the regularity condition specifies w'(0) = w"(0) = (see Ea. (|?T|l ). 

Unfortunately, under these boundary conditions we could not find soliton solu- 
tions in Eq. (|^ for any value of 52 ■ 

Recently, Forkel studied the soft mode action of the Yang-Mills the ory and 
investigated the solitonic, saddle point solution with hegdehog symmetry His 
work is quite encouraging to find the class of solutions in the infrared part of the 
Yang-Mills action. But, the solutions the author obtained are solitons with non- 
zero topological charge and, in the case of the Hopf solitons which possess zero 
topological charge, the situation is a little more complicated. 

From the identity 

j d^x[{d^n ■ d^n) - [df^nY] ^ J d^x{d^n x n)^ , (24) 
one easily finds that the static energy obtained from the last two terms in Ea. pOfl 

) = J d?x{d^n X nf (25) 

gives the positive contribution. The total static energy is stationary at = 

\J E^/iE^l'' + E^^'') and hence the possibility of existence of soliton solutions is 
not excluded. And also, the positivity of Ea. (|25|l does not spoil the lower bound Q 
of original SEN and the possibility still remains. 

Therefore, we suspect that the failure of finding the stable soliton is caused by 
the fact that higher derivative theory has no lower hound. We shall investigate the 
lower bound in the higher derivative theory in detail in the next section. 

5. Higher derivative theory 



In this section, we address the basic problems in the higher derivative the- 
oryEl 22 23 24 25 -^yj^jf.}^ essentially falls into two categories. The first problem con- 
cerns the increase in the number of degrees of freedom. Eor example, if the theory 
contains second derivative terms, the equation of motion becomes up to the order in 
the fourth derivative. Thus, four parameters are required for the initial conditions. 
If one considers more higher order terms, the situation gets worse. However, this is 
not serious problem for our study because our concern is the existence of static soli- 
ton solutions. The second problem is that the actions of the theory are not bounded 
from below. This feature makes the higher derivative theories unstable. 

The lagrangian and the hamiltonian formalism with higher derivative was firstly 
developed by Ostrogradski We consider the lagrangian containing up to nth 
order derivatives 

I dtC{q,q,--- ,q^"^). (26) 
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Taking the variation of the action SS = leads to the Euler-lagrange equation of 
motion 

i:<-')'f(iT)=°- ™ 

i=0 ^ 



The hamihonian is obtained by introducing n generahzed momenta 



E(-l)^-^£r^(7§T)' - V-,n, (28) 

0=1+1 

or 

dC dC d 

P-^d^^ P^^d^-d-tP^^'^ z = (29) 

and n independent variables 

qi = q, q, = q^'-^\ i = 2,--- ,n. (30) 

The lagrangian now depends on the n coordinates qi and on the first derivative 
(fn — q^^^^ ■ The hamiltonian is defined as 

n n — 1 

'H{qt,Pi) = ^Piqi - C = ^Ptqt+i + PnQn - >C . (31) 

i=l i=l 

The canonical equations of motion turn out to be 

9^ . dn 
Qi^T^, Pt = --^- 32 
opi oqi 

Thus, we replace a theory of one coordinate q system obeying 2n— th differential 
equation with a set of 1-st order canonical equations for 2n phase-space variables 

We consider a simple example including a second derivative terml^, defined as 

£ = i(l + eW-^-V-^eV, (33) 

where constant e works as a coupling constant of second derivative term. The equa- 
tion of motion is 

{l + e^Lj'^)q + uj^q + e^q'^'^'^ ^Q. (34) 

From Eq. if^ . one gets 

2- dC d fdC\ 2 2\ ■ 2--- 



dq ' dq dt\dqJ 

Thus the hamiltonian becomes 



H — XTTq + ij-Kg — C 
2e 



- T^^^l - o (1 + + 0^^^^ ■ (36) 
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We introduce the new canonical variables 



UJ 



1 / 2 2 • \ '"^ / \ 

-.(e UJ q-TTg), p+= =={q-TTq). 



vl — £Vl — 



and the hamiltonian can be written in terms of these variables as 

The corresponding energy spectra is then given by 

E={n+l)LO-im+^)- , 7i,m = 0,1,2,- •• (37) 
2 2 £ 

One can see that in the limit e — > the energy goes to negative infinity rather than 
approaching to the harmonic oscillator energy eigenstates. 

To obtain physically meaningful solutions, we employ the perturbative analysis 
where the solution is expanded in terms of the small coupling constant and the 
Euler-Lagrange equation of motion is replaced with the corresponding perturbative 
equation. The solutions of the equations of motion that are ill behaved in the limit 
e — *■ are excluded from the very beginning . 

We assume that the solution of Eq. (|34|l can be written as 

oo 

9pcrtW-5]e"gW. (38) 

n=0 

Substituting Ea. (|38|l into Ea. (|34|l and taking time derivatives of these equations, 
we obtain the constraints for higher derivative terms 

0(e") 

equation : % + t^'^'Zo — , (39) 
constraints : Qo = ^^'^'io, 9 o = '-^'^Qo ■ (40) 

equation : q2 + io'^qo + uo'^q2 + 9 o = , 

^ q2 + uj'^q2 = , (using (jUJ) , (41) 
constraints : 9 2 = —^JJ^q2^ i 2 — (^*q2 ■ (42) 

equation : q4 + uj'^q2 + uj^q^ + 9 2 = 0, 

^ qi + u'^qi = , (using (O) , (43) 
constraints : q 4 — —0^^(74, 9 4 — ^^^94 • (44) 
Combining these results, we find the perturbative equation of motion up to O(e^) 

gport + w^gpcrt = O(e^) • (45) 
which is the equation for harmonic oscillator. 
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6. Search for the sohton solutions (2) perturbative analysis — 

Let us employ the perturbative method introduced in the last section to our prob- 
lem. We assume that 52 is relatively small and can be considered as a perturbative 
coupling constant. Thus, the perturbative solution is written by a power series in 

92 

00 

w{7j)^J29>niv)- (46) 

ra=0 

Substituting Ea. (|46|) into Ea. H23|l . we obtain the classical field equation in 0(172) 

fo{wo,WQ,Wg) + gifi{wo,w'Q,WQ) = 0. (47) 

Taking derivatives for both sides in Eq. (|47l) and solving for Wq , wl^"^ , w^^"' read the 
following form of the constraint equations for higher derivatives 

u;« =i^W(«,o,u.^), i = 2,3,4. (48) 

The equation in 0(52) can be written as 

{fo+92fi)oigl)+f2{w,,wiw'^,wi'\wi^^)=0. (49) 

Substituting the constraint equations H48|l into Ea. (|49|l and eliminate the higher 
derivative terms, one can obtain the perturbative equation of motion 

fo{w, w', w") + gihiw, w',w") + .92/2 (w, w') = 0(gl) . (50) 

One can sec that Ea. (|50() has stable soliton solutions. 

The numerical results of the hmctions w{rj) are displayed in Fig[21 In Fig|21 we 
show the energy density plot for the original SEN model in the cylindrical coordi- 
nates {p, z). In Fig^ we present the energy density for our extended soliton in the 
(p, z). Both results share the toroidal shape and no notable difference at least for the 
small couphng constant g2- In Fig|Sl we plot the total energy for H = l,2,gi = 0.4 
as a function of the coupling constant g2- As can be seen, the change is moderate 
and, more interestingly, the energy seems to be linearly dependent on 52 at the 
region of smaller 52- For larger 52, the data gradually deviate from the linear be- 
havior. The deviation would be due to the fact that our analysis is based on the 
first order perturbation. Therefore, it is possible to observe the critical values of g2 
for each gi and H in which our simple first order perturbation is valid. Performing 
linear fitting, one finds 



Eh=i - 131 + 121(72 , 
Eh=2 208 + 102^2 . 
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Fig. 2. The function w{x) of (a) H = l,g]_= 0.4, 32 = 0.25, (b) H = 2,gi = 0.4,^2 = 0.3 (the 
rescahng radial coordinate x = rt/(l — rf) is used), together with the results of the original SFN 
model. 
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Fig. 3. The energy density of naive SFN model in the cyhndrical coordinates (p, z) for H - 
2,gi = 0.4. 



In Figl^l we plot the energies for various values of gi as a function of g2- They are 
fitted to the following linear functions 

^51=0.2 - 147 + 14252 , 



91=0.6 
'ff 1=0.8 



255 + 8352 : 
294 + 7252 . 



From these results we can extract formulation of topological bound for the energy 
of the second derivative term -E2nd: 



V9i 



(51) 



If /3 ~ 76 is chosen, all data are well fitted by this single formula. For the energy 
from standard SFN action, we tentatively employ the topological lower bound l|3Jl. 
Then the topological lower bound including second derivative terms becomes 



(52) 



with a — 167r^. Though this formula is not based on any theory, it contains in- 
teresting information about the spectra of soliton. That is , l|52(l has the minima 
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for given (72- The energy bound is thus Ogtt — "^cty/^H^/^ . 

Let us summarize the numerical results presented in this section. If we take into 
account the second derivative term, the obtained soliton mass will become twice 
value of the one in the naive SFN soliton. In FigEl a schematic plot of Ea. H52|l is 
shown, and the existence of such local minima is observed. Interestingly, though 
at first glance our lower bound formula 1)52(1 has no 3/4-scaling behavior like 
essential to knotted solitons the energy at local minima recovers this scaling 
law. 

Unfortunately, such minima can not be observed at present within our numerical 
framework. Because, from the condition ((53|l . the value of (72 should be almost twice 
of gi for iJ = 1, and four times for H — 2. Otherwise, we need to take into account 
the higher order perturbation terms to employ totally different formalism to achieve 
these values. 
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Fig. 5. The static energies of the sohton for H = 1, 2 as a function of g2 {gi = 0.4). Also 
the Unearly fitted Une of the H = 1 : Egtt = 131 + 121^2 (dotted hne) and of the H = 2 : 
Bstt = 208 + 102^2 (dashed Une). 



7. Summary 

In this paper we have studied the Skyrnie-Faddeev-Niemi model and its extensions 
by introducing the reduction scheme of the SU(2) Yang-Mills theory to the corre- 
sponding low-energy effective model. To simplify the matter, we proposed an ansatz 
for n. That extensively reduced the computational time and did not affect to the 
property of the soliton solution. The requirement of consistency between the low- 
energy effective actions of the YM and the SFN type model leads us to take into 
account second derivative terms in the action. However, we found that such an ac- 
tion including the second derivative terms does not have stable soliton solutions. 
This is due to the absence of the energy bound in higher derivative theories. This 
fact inspired us to employ the perturbative analysis to our effective action. Within 
the perturbative analysis, we were able to obtain soliton solutions. 

It should be noted that our solutions do not much differ from the solution 
of original SFN model, at least in the perturbative regime. We suspect that an 
appropriate truncation (for instance "extra fourth order term + second derivative 
term") can supply stable solutions that are close to the original SFN model. Thus 
we conclude that the topological soliton model consisting of the "kinetic term + a 
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Fig. 6. The static energies of the soliton for gi = 0.2,0.4,0.6,0.8 as a function of g2 (H = 2). 
Also the Hnearly fitted Une of the gi = 0.2 : Estt = 147 + 142^2 (dot-dashed hne), the gi = 0.4 : 
Estt = 208 + 10232 (dotted line), the gi = 0.6 : Estt = 255 + 83g2 (dashed line) and of the gi = 0.8 
: Sstt = 294 + 7292 (solid line). 

special fourth order term" like SFN model is a good approximation. 

Our analysis is based on perturbation and the coupling constant ^2 is assumed 
to be small. However, Wilsonian renormalization analysis of YM theory suggests 
that the coupling constants 51,52 (and the mass scale parameter A) depend on 
the renormalization group time t = logfc/A (fc,A are infrared, ultraviolet cutoff 
parameter) and they are almost comparable. To improve the analysis, we could 
perform the next order of perturbation, but it is tedious and spoils the simplicity 
of the SFN model unfortunately. 

We found numerically that the energy of the extended soliton solutions is linear 
in the coupling constant 52 and then extracted a new mass formula for the soliton 
solutions including second derivative term. We expect that the global minima for the 
coupling constant 52 exists and the corresponding energy becomes twice the value of 
the naive SFN action. Of course, this statement is not based on any theory and we 
have to wait for its theoretical confirmation. Since our mass formula was obtained 
from numerical study in the perturbative approach, it is uncertain whether such 
linear behavior is kept for larger coupling constant 52 {like twice of gi for H = 1 
and four times for H = 2). To confirm that, we should proceed to investigate next 
order perturbation, or, otherwise, find some analytical evidence of this formula. We 
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Fig. 7. The plot of the conjecture 1521 for (gi, 92)1 in the case of = 1. 

point out that the perturbative treatment is only used for excluding the ill behavior 
of the second derivative field theory. We believe that applying this prescription 
should not alter the essential feature of the soliton solutions, e.g., existence of the 
solutions, linear behavior of the mass spectra, et.al.. 

Finally, let us mention the application of the soliton solutions to the glueball. 
Obviously this is one of the main interest to study the model and, many authors 
have given discussions about this J^^^o^-^. On the other hand, the possibility of 
the magnetic condensation of the QCD vacuum within the Cho-Faddeev-Niemi- 
Shabanov decomposed Yang-Mills theory have been studied by Kondcl^. The au- 
thor claims the existence of the nonzero off diagonal gluon mass Mx, which is 
induced in terms of the condensation of the magnetic potential part of the decom- 
position ^ {dfiTi) X n, as 

Afi = {B^ ■ B^) = {{d^nf) . (54) 

Throughout our calculation, we set the overall coupling constant A = 1 in the action 
(|T]l but, in this sense, it should reflect the information of such gluon mass, or the 
condensation property of the vacuum. After a careful examination of the value of 
the property of A, we will be able to accomplish the detailed predictions for the 
glueball mass. 
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